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Abstract 

A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are 
developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. 
This surrogate fuel model depicts a fuel blend that properly represents the general physical and 
chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the 
surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series 
of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model 
consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation, 
and a detailed wet-CO mechanism. The final thermophysics formulations are incorporated with a 
computational fluid dynamics model for prediction of the combustor efficiency of an uni-element, 
tri-propellant combustor and the radiation of a kerosene-fueled thruster plume. The model 
predictions agreed reasonably well with those of the tests. 

Nomenclature 

A = pre-exponential factor 
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Cp = heat capacity 

E = activation energy 

H = enthalpy 

K = forward rate constant 

P = pressure, atm 

Q = radiative heat flux, , btu/ft -s 
R = universal gas constant 

S = entropy 

T = temperature, K. 

a = species mass fraction 

0 = view angle, deg 

Subscripts 

c = combustion or chamber 

f = formation or forward reaction 


Introduction 

Russian-built kerosene-fueled rocket engines such as RD-170* or its U.S. proposed 
counterparts such as RD-704 and NASA MSEC developed Fastrac engines have been identified as 
potential candidates to fly the Single-Stage-to-orbit Reusable Launch Vehicles. In order to support 
the associated engineering issues, specifically the preliminary conceptual design and evaluation of 
the performance of the injectors and thrust chambers using Computational Fluid Dynamics (CFD), 
accurate and computationally efficient models that properly represent the thermophysics, namely 
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fuel formula, thermodynamics, and finite-rate combustion-kinetics, have to be used. Unfortunately, 
models pertaining to these aspects were underdeveloped. 

In this study, based on reported physical-chemical property data, a one-formula surrogate fuel is 
proposed as a generic representation for kerosene, or its derived fuel RP-1. The thermodynamics of 
the surrogate fuel such as heat capacity, enthalpy and entropy are generated and the rationale for 
generation is given. A kerosene combustion kinetics is proposed based on a quasi-global kinetics 
format and the rationale for which, including the specific reactions chosen, is also given. The final 
thermophysics characterization is implemented in a CFD model and benchmarked on two rocket 
applications: the combustion efficiency of an experimental tri-propellant combustor and the 
radiation of a test rocket plume. 

One-Formula Surro2ate Fuel Model 

Common hydrocarbon rocket fuels such as kerosene, or jet fuels are derived from petroleum, 
while RP-1 is a straight run kerosene fraction^’^. These hydrocarbon fuels are complex mixtures 
of many components and their exact composition and properties vary from batch to batch. Some 
information on these substances has been reported in the literature. For example, it is reported'* 
that kerosene typically consists of ten hydrocarbons containing 10-16 carbon atoms per 
molecule, while the constituents include n-dodecane, alkyl derivatives of benzene, naphthalene 
and its derivatives. Also, at least 87 identifiable hydrocarbons are reported in RP-1^. 

It is obvious that using all 87 identifiable hydrocarbons to represent RP-1, or even ten species 
for kerosene would be highly inefficient in any kind of performance calculations. A surrogate 
fuel model composed of a neat compound or several neat components for kerosene-derived fuels 
has to be developed, while important thermo-physical properties are retained. Such a 
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simplification is necessary for computation-intensive design calculation using complex fuels or 
blends. Conventionally, elemental formulas such as CHi.9423 and its derived standard heat of 
formation have been used in one-dimensional theoretical rocket performance calculations*’ for 
JP- 4 /RP -1 powered engines. However, elemental formulas have severe limitations when design 
calculations other than one-dimensional performance is desired, especially those using CFD 
codes where molecular fuel formulas are preferred. In the past, several non-elemental one- 
formula surrogate fuel models were proposed. For example, one-formula models CioHi9^ and 
Ci 2H23‘^ were used to represent kerosene, whereas Ci2H26*’^ was used to replace RP- 1 . The 
simplistic nature of these one-formula fuel models makes them easy to use. However, these 
simple formulas also present subtle problems when used in performance calculations. For 
example, when comparing to the physical-chemical properties of kerosene and RP- 1 , the 
molecular weight of CioHjg is too low, the H/C ratio of C12H23 is too low and that of C12H26 is 
too high. Furthermore, C12H26 (n-dodecane) is a paraffin while kerosene normally contains only 
41 % non-cyclic hydrocarbons. 

Other than the elemental formula and one-formula fuel models, a more detailed modeling 
approach is the multiple-formula models. For example. Farmer and Anderson^ used a three- 
formula surrogate fuel model to represent RP- 1 , i.e., 17 . 4 % C13H12 (methylbiphenyl), 45 . 4 % 
C12H24 (n-heptylcyclopentane), and 37 . 2 % C12H28 (n-tridecane). With multiple molecules, the 
physical-chemical properties of kerosene or RP -1 may be matched better. However, this 
approach often creates extra number of chemical species thereby slowing down the computation. 
Weighing the advantages and disadvantages of the one-formula and multiple-formula 
approaches, a one-formula surrogate fuel C12H24 is hereby proposed as a generic representation 
for kerosene or RP- 1 , mainly for computational efficiency reason since only one parent fuel is 
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involved. In addition, although C 12 H 24 takes the form of a one-formula surrogate fuel, it is 
modeled (kinetically) as 41.7% paraffin and 58.7% naphthene, thereby matching the approximate 
paraffin/naphthene split of 41/56.^’^ In fact, C 12 H 24 does not represent a certain molecule, rather 
a mixture of many neat components that has the same averaged thermo-physical-chemical 
characterization as that of the reported kerosene/RP-1. Therefore C 12 H 24 has the advantages of a 
multiple-formula surrogate fuel without being represented as multiple fuels. In summary, the 
(averaged) molecular weight, elemental formula and formula weight, and the paraffin/naphthene 
split of the surrogate fuel C 12 H 24 match reasonably well with those of reported kerosene and RP- 
1, as shown in Table 1. The paraffin-to-naphthene split is an important factor to match in terms 
of accurate soot prediction since soot forms much easier from naphthene than from paraffin. 
This point will be elaborated in the quasi-global combustion kinetics section. Notice the small 
amount of aromatics and olefins in kerosene and RP-I are ignored and lumped into the paraffins 
and naphthenes of the surrogate fuel model. 

Table I. Comparison of thermo-physical-chemical characterization of model fuel with reported 
data. 



Kerosene 

RP-1 

Surrogate fuel 

Molecular Formula 



C12H24 

Molecular Weight 

175'° 

172-175^ 

168 

Elemental Formula 


CH,95-CH2.o", 

CH20 



^TT 6 

1.9423 1.953 


Formula Weight 


13.97-14.03^13.97° 

14.03 


5 


He, kcal/g 

-10.278\ 

- 10.24 1\ 

-10.278 


-10.321^ 

-10.356^ 


Hf, 298 K. kcal/mole 


-5.430‘^/CH, 9423 

-92.200/C, 2 H 24 , 




- 7 . 683 /CH 20 

Cp,5i6k,i atm cal/molc-K 


101^ 

103 

Paraffins (n & iso) % 


41^ 

41.7 

Naphthenes % 


56^ 

58.3 

Aromatics % 

5^ 

5",3^ 

0 

Olefins % 

1 ' 

0^ 

0 


Gaseous-Phase Surrogate Fuel Thermodynamics Generation 

An important thermodynamic property to be considered for the surrogate fuel is the heat of 
formation if accurate heat release during combustion is to be predicted. In general, the heat of 
formation can not be measured directly but must be determined indirectly from its heat of 
combustion, measurable calorimetrically. The heat of formation for the surrogate fuel is 
therefore determined by matching its heat of combustion with that of kerosene/RP-1. This is 
accomplished by writing a complete combustion reaction for C12H24: 

C 12 H 24 + 18 O 2 ^ 12 CO 2 + 12 H 2 O 
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The heat of reaction of the above equation is the heat of combustion. By definition, the net heat 
of combustion is obtained when the product H 2 O is at its gaseous state. The heat of formation of 
the surrogate fuel is written as 


Hf,C12H24 =12 Hf,C 02 + 12 Hf ,H20 - He 

The above equation underscores the importance of an appropriate fuel model since the 
stoichiometry of the combustion reaction determines the value of the molar heat of formation and 
eventually the amount of combustion products in the engine system. The simplicity of the one- 
formula surrogate fuel model makes it unnecessary to deal with multiple surrogate fuels, but 
directly from the reported heating value of kerosene/RP-1. Several reported standard heats of 
combustion are listed in Table 1. The maximum difference among them is only 0.115 kcal/g and 
amounts to about 1 % of the heating value. These values are deemed as consistent and are within 
the uncertainty bound of the measurement. A midrange value of -10.278 kcal/g is chosen. As a 
result, a thermodynamically consistent heat of formation of -92.200 kcal/mole is derived for 
C 12 H 24 and a value of -7.683 kcal/mole is calculated for its elemental form CH 2 0 i shown in 
Table 1. 

The next step is to construct the three thermodynamic functions of heat capacity, enthalpy, 
and entropy as functions of temperature in a usable form. The standard fourth-order polynomial 
heat capacity form*^ is used: 
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R 


-^ = a, + 




RT 


2 3 

“ 2 ^ 


-J- = a.\nT-¥a^T + ^T^ 

R ^ 2 2 3 4 


r 


+ Cl'j 


The enthalpy of the surrogate fuel is constructed using the heat capacities of n-deodecane, due to 
the closeness of n-dodecane (C 12 H 26 ) with C 12 H 24 family of molecules in terms of carbon and 
hydrogen atomic numbers and the notion"*’^ that n-dodecane is a major component in 
kerosene/RP-1. Notice the enthalpy (heat) of formation is anchored with the thermodynamically 
consistent heat of combustion of kerosene/RP-1. Hence, the heating curve generated herein is 
more realistic than that of Ref. 4 in which the enthalpy of n-dodecane was used as a placeholder. 
The generated heat capacity for the surrogate fuel matches well with that of RP-1, as indicated in 
Table 1. The entropy of formation of n-dodecane, is used to construct the entropy curve for 
C 12 H 24 . The resulting least square coefficients fitted for two temperature ranges are listed in 
Table 2. 


Table 2. Thermodynamic coefficients for C 12 H 24 . 



1000 to 5000 K 

300 to 1000 K 

ai 

0.36440206E+02 

0.39508691E-1-01 

a2 

0.54614801E-01 

0.10207987E-I-00 

as 

-0.16091 151E-04 

0.13124466E-04 

a4 

0.21478497E-08 

-0.76649284E-07 

as 

-0.10131 180E-12 

0.34503763E-10 
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-0.63890109E+05 


-0.52093574E+05 


a6 


a? 


-0.15798973E+03 


0.2198095 lE+02 


Notice that the reported magnitude of the heat of vaporization of kerosene (0.059 kcal/g at 
normal boiling point)^ is not only much smaller than that of the heat of combustion, but is also 
less than the error bound (0.115 kcal/g)^’^’^’^. In addition, the latent heat decreases to zero at 
critical pressure, as most of the rocket engines are operated at higher pressures. These rationales 
allow the gaseous-phase thermodynamics to be curve-fitted to 300 K, about 142 K lower than the 
normal boiling point (542 K)^ of kerosene. Nevertheless, in actual calculations involving liquid 
kerosene fuel, the heat capacity and latent heat of vaporization of liquid kerosene‘s are used for 
the liquid-phase thermodynamics, as shown in latter sections of this study. 


Surrogate Fuel Thermodynamics Validation 

From a thermochemical consistency standpoint, the proposed one-formula surrogate fuel and 
its thermodynamics can be validated by performing a series of one-dimensional theoretical 
rocket engine performance calculations, using the thrust chamber specifications and operating 
conditions of a Russian engine RD-no’ '®. Fig. 1 shows a comparison of the calculated 
theoretical chamber and nozzle exit temperatures as a function of mixture ratios for C 12 H 24 and 
its elemental formula CH 2.000 for the self-consistency test. The calculated temperatures for 
C 12 H 24 and CH 2.000 coincide, as expected, since their heat of formations are anchored with the 
chosen heat of combustion at -10.278 kcal/g. Also shown in Fig. 1, are the results using 
elemental formula CH] 9423 with heats of formation of -5.430, -6.613, and -7.133 kcal/mole. 
These values correspond to heats of combustion of -10.356^, -10.278^, and -10.241^ kcal/g. 
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respectively. The difference among the curves is small since the difference in heats of 
combustion is small, especially in the fuel lean region. When the temperature profiles are 
compared at the same heat of combustion (- 10.278 kcal/g), the minor difference between CH2 
and CHi. 9423 curves is the result of their difference in H/C ratios. 

Fig. 2 shows the predicted theoretical nozzle exit gas composition. Again, the mole fractions 
resulted from using C12H24 coincide with those of CH2.ooo- Graphite carbon (Cs) and methane gas 
formed at nozzle exit at low mixture ratios due to fuel rich combustion at low exit temperatures, 
as expected. These results show that the proposed one-formula surrogate fuel model is suitable 
for one-dimensional theoretical rocket performance calculations and is thermodynamically 
consistent. However, it should be noted that the system-engineering type of calculations often 
overpredict the performance, especially for heavy hydrocarbon fueled engines. This is because 
performance impacting thermophysical processes such as ignition delay times and vaporization 
are usually not considered. In addition, other performance impacting and design particular 
physical processes such as mixing, shock losses, geometry losses, film cooling and boundary 
layer losses are also not included. And those are motivations for this study. 

Ouasi-Global Finite-Rate Combustion Kinetics 
Kerosene Decomposition and Wet-CO Mechanism 

Detailed kinetics mechanism involving elementary steps for parent fuel decomposition is 
computationally prohibitive for considering kerosene/RP -1 combustion in CFD calculations. On 
the other hand, one-step'* or multiple-step global kinetics models'^ is computationally efficient, but 
frequently over-predicts the flame temperature, or not generic enough to take advantage of the 
modem diagnostic techniques. For example, hydroxyl radical (OH) is usually not involved in those 
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models and laser-induced fluorescence imaging of the OH-radical flame structure data"^ can not be 
utilized. The quasi-global kinetics that combines several global steps with a detailed wet-CO 
mechanism is probably best suited for describing complex fuel combustion in a computationally 
intensive environment. The original quasi-global kinetics*^ was established based on the 
observation that straight-chain hydrocarbons and cyclic hydrocarbons have distinctively different 
ignition delay time characteristics, and that the ignition delay times for straight chain 
hydrocarbons are similar and so are those for the cyclic hydrocarbons. That important 
observation led to the development of the original quasi-global kinetics model in which two 
irreversible global steps are used to describe the decomposition of the straight chain 
hydrocarbons and the cyclic hydrocarbon into intermediate products, whereas a detailed wet-CO 
mechanism converts the intermediates into the final combustion products. In the quasi-global 
kinetics format, OH radical is not only considered (in the wet mechanism), but plays a vital role 
in the accurate prediction of the heat release during hydrocarbon combustion. The original quasi- 
global kinetics scheme'^ was evaluated with a variety of simplified flow configurations, 
including the stirred reactors, plug flow reactors, and turbulent diffusion flames. It was later 
expanded to an “extended” quasi-global kinetics model in which two more intermediate 
species and several more global kinetic steps are added to better describe the fuel-rich 
combustion. In this study, under the premises of computational efficiency and the framework of 
one-formula surrogate fuel model, the original quasi-global kinetics format is considered. 
Following that format, two global steps are proposed: one for the paraffin portion and another for 
the naphthene part of the surrogate fuel. The rates of the two global steps are modified directly 
from those of the straight chain and cyclic global steps'^ according to the paraffin and naphthene 
split in the proposed surrogate fuel (C12H24), thereby satisfying the ignition delay times 
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requirement of the original quasi-global kinetics model. The proposed kerosene global steps and 
the existing CO-wet mechanism are shown in Table 3. The “standard form” in the fourth column 
of Table 3 means the species concentration terms in the rate expression follows the stoichiometry 
of the reaction. Also shown in Table 3, the forward reaction rate constant is expressed in the 
standard Arrhenius form. The backward reaction rate of the reversible reaction is calculated 
from its forward rate and equilibrium constant. 

Soot Formation 

Under fuel rich conditions, kerosene/RP-1 forms soot readily. This is because Naphthene and 
aromatic hydrocarbons form soot rapidly (condenstaion-polymerization) by directly condensing 
themselves into polycyclic aromatic hydrocarbons (PAH). On the other hand, paraffins form soot 
slowly. This is because paraffins have to be break up into smaller fragments first, from which 
fusing of the fragments occurs to form naphthenes and aromatics, and PAH’s form eventually and 

indirectly (fragmentation-polymerization).’^ '^ These PAH’s, also known as soot precursors, are 

20 

then undergone a series of physical processes to form coagulated soot particles. Frenklach et al. 
developed a comprehensive soot formation mechanism in which 180 species and 619 elementary 
reactions are used in an attempt to describe the aforementioned soot formation processes. 
However, at the present moment, it is far too expensive to be incorporated into a CFD code while 
the oxidation of those 180 species was not even considered. Nickerson and Johnson^ abridged 
Frenklach et al.’s model to 19 reactions for prediction of the soot formation in LOX/RP-1 injectors. 
However, their model completely ignored the fast condensation-polymerization process, and a soot 
formation mechanism of 19 reactions is still of considerable size. A global step that makes soot 
directly from C 12 H 24 is used instead, as shown in Table 3. For convenience, pseudo-gas graphite 


12 



carbon (Cs) is used to represent soot. The forward reaction rate is modified from that established 
for toluene combustion in jet-stirred combustors and shock tubes. The effects of both 
condensation-polymerization and fragmentation-polymerization are included in this global step. 

Heterogeneous Soot Oxidation 

The rates of carbon oxidation were first measured and fitted with a heterogeneous reaction 
model involving three elementary steps and two active carbon sites to forming CO. These rates 
were adopted by several investigators as the rates for heterogeneous soot oxidation. For 
example, earlier studies'*’’^* used a one-step equation that reacts soot with molecular oxygen to 
form CO 2 , whereas the reaction rates are those of Ref. 22 but the rate expression was modified 
slightly. In 1995, Hier^^ reported that Roth et al.^"* and Olander et al.^^ independently confirmed 
that the reaction forming CO is much more likely than that forming CO 2 . Roth et al.^“* used laser 
absorption to confirm the presence of CO and the absence of CO 2 , whereas Olander et al. used 
mass spectrometry to determine that the reaction forming CO was at least two orders of 
magnitude more probable than that forming CO 2 . Based on those assertions, this study simplifies 
the three-step reactions^^ to a one-step global kinetics forming CO, while maintaining the 
original rates and the heterogeneous rate expression form^^, as shown in Table 3. 

Homogeneous Soot Oxidation 

Not only is the OH radical an important controlling species of the heat release in combustion 
processes, as noted by one of the major premise of the quasi-global kinetics. There is evidence that 
OH dependent oxidation of the soot needs to be considered under many flame conditions, 
especially in hydrocarbon fueled exhaust plumes.^*^'^^ The OH dependent homogeneous soot 
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oxidation reaction of Slack et al.^^ is included in the kerosene quasi-global combustion kinetics so 
that the effect of soot thermal non-equilibrium on plume characteristics can be assessed, as shown 
in Table 3. Notice this reaction is endothermic, whereas the heterogeneous soot oxidation reaction 
is exothermic. Also, in stead of the irreversible reaction in its original form,^*^ a reversible reaction 
is used in this study since it produces better results with the following validation cases. 


B -E/RX 

Table 3. Kerosene/RP-1 quasi-global combustion kinetics mechanism. Kf = AT e' 


Reaction 

A 

B 

E/R 

Form 

Ref. 

Paraffin Global Step 

C 12 H 24 + 6 O 2 ^ 12CO - 1 - I 2 H 2 

3.888E4 

1 

1.220E4 

p"'[C,2H24]”'[02] 

This Work 

Naphthene Global Step 

C 12 H 24 + 6 O 2 — > 12CO + I 2 H 2 

2.312E7 

1 

1.965E4 

p"'[C,2H24]"'[02] 

This Work 

Wet-CO Mechanism 

H 2 + O 2 — OH + OH 

1.700E13 

0 

2.407E4 

Standard 

15,16 

OH - 1 - H 2 = H 2 O + H 

2.190E13 

0 

2.590E3 

Standard 

15,16 

OH + OH = 0 - 1 - H 2 O 

6.023E12 

0 

5.500E2 

Standard 

15,16 

0 - 1 - H 2 = H + OH 

1.800E10 

1.0 

4.480E3 

Standard 

15,16 

H -h O 2 = 0 - 1 - OH 

1.220E17 

-0.91 

8.369E3 

Standard 

15,16 

M-f-0-hH = OH + M 

1.000E16 

0 

0 

Standard 

15,16 

M + 0 + 0 = O 2 + 

2.550E18 

- 1.0 

5.939E4 

Standard 

15,16 


14 


M + H + H = H2 + M 

5.000E15 

0 

0 

Standard 

15,16 

M + H + OH = H 2 O + M 

8.400E21 

-2.0 

0 

Standard 

15,16 

CO + OH = H + CO 2 

4.000E12 

0 

4.030E3 

Standard 

15,16 

CO + O 2 = CO 2 + 0 

3.000E12 

0 

2.500E4 

Standard 

15,16 

CO + 0 + M = CO 2 + M 

6.000E13 

0 

0 

Standard 

15,16 

Soot Formation Global Step 






C,2H24^ 12Cs+ 12Hz 

4.4947E15 

-1.94 

1.610E4 

[C,2H24]‘''[02]°' 

This Work 

Heterogeneous Soot Oxidation 






Cs + O. 5 O 2 ^ CO 

1.0 

0. 

0. 

72Rox[Cs]/(PsDs) 

This Work 

Rox = KaP02X/(1+KzP02) 





22 

+ KbPo2(1-X) 






%= 1/(1+(Kt/Kb)Po2) 






Ka 

2.0000E1 

0 

1.5098E4 



Kb 

4.4600E-3 

0 

7.6497E3 



Kt 

1.5100E5 

0 

4.8817E4 



Kz 

2.1300E1 

0 

-2.063E3 



Homogeneous Soot Oxidation 






Cs + OH = CO + H 

1.2200E9 

0.5 

0 

Standard 

This work 
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* M stands for third-body collision partner 


Application to Kerosene fueled Thrust Chambers 
Pennsylvania State University Tri-Propellant Combustor 

The maximum combustion chamber pressure achievable is usually associated with the 
combustion efficiency. The proposed thermophysics model is incorporated with a CFD code, 
FDNS^^’^^, to compute the chamber pressure of an uni-element shear tri-axial tri-propellant 
combustor. The test rig comprises of several sections, including an injector assembly, igniter, 
window and blank sections, and a nozzle assembly. Liquid RP-1 (LRP) is introduced through 
the central post of the tri-axial tri-propellant injector and swirled using a tangential swirl nut, 
whereas gaseous H 2 (GH 2 ) and gaseous O 2 (GO 2 ) are injected through the first and second 
(outer) annuluses, respectively. Details of the injector and the thrust chamber setup can be found 
in Refs 14 and 30, whereas the chamber operating conditions used for model validation can be 
found in Ref. 30. 

FDNS is a multidimensional, multispecies, viscous flow, pressure-based reacting flow solver. 
It was developed at Marshall Space Flight Center (MSFC) and is continuously being improved 
by MSFC personnel and its supporting contractors. By incorporating and validating the 
proposed thermophysics formulation of kerosene combustion with a flow solver like FDNS, 
realistic performance computations can be performed to support the design of kerosene fueled 
injectors and thrust chambers, rather than relying on one-dimensional type of system 
performance calculations. FDNS solves simultaneous liquid-droplets-gas dynamics by 
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combining the volume-of-fluid (VOF) and Eulerian/Lagrangian tracking methods into an unified 
algorithm for efficient calculation of multiphase flows at all speeds. The gas-liquid interface 
mass, momentum and energy conservation properties are modeled by continuum surface 
mechanisms. Details of the solution procedure for liquid-droplets-gas governing equations can be 
found in Ref. 29 and 31. In the VOF formulation, the weight-averaged conservation equations of 
mass, momentum, and scalar variables are formulated in an Eulerian framework. Current 
formulation considers the vaporization of liquid droplets into gases, but not the direct 
vaporization of liquid jet core into gases. An empirical mass stripping rate is applied to the VOF 
equation along the liquid-gas interface for liquid intact core resolution. The stripped liquid mass 
undergoes secondary breakup process to form smaller droplets. Lagrangian particle tracking 
method is used to treat the dynamics and heat/mass transfer of droplets, as computational parcels 
in a statistical sense. Energy and momentum are solved for each parcel for every species and 
relative enthalpy is the dependent variable for the droplet energy equation. For simplicity, parcel 
temperature is assumed to be constant at any time instant, but varies throughout its life 
expectancy via local heat transfer at the surface and vaporization. The heat capacity and latent 
heat of vaporization for liquid droplets are taken from Ref. 12. 

An axisymmetric thrust chamber is formulated and the computational domain covers the 
entire thrust chamber. As such, there is no ambiguity in the exit flow boundary condition once 
the throat flow is choked. For comparison purpose, a gaseous RP-1(GRP)/GH2/G02 injector 
flowfield is computed first. After its convergence, the central-post flow passage is then replaced 
with LRP and a LRP/GH 2 /GO 2 injector flowfield computation is ensued. To start the calculation, 
a cold-flow thruster flowfield is computed first until the nozzle flow at the end of the combustor 
choked. The propellants in the mixing layers are then ignited by temporarily imposing hot spot 


17 



with an elevated temperature (1000 K) in the flowfield. Once ignited, the flame propagates 
upstream and spreads downstream, whereas the final solution is not influenced by the location 
nor the added energy of the temporary hot spot. Due to the backward-facing step formation 
created by the injector faceplate and the side-wall, a flow recirculation forms to serve as a flame 
holder. A fixed mass flow rate boundary condition is imposed at the flow inlet and the inlet total 
enthalpy is conserved. The flow properties at the nozzle exit plane are extrapolated from those of 
the interior. 

Figure 3 shows the computed scalar contours for the unielement GRP/GH2/GO2 
injector/combustor. Since the LRP Jet is replaced with the GRP jet, there is no oClrp contours. The 
first three scalar contours are those of the reactants GRP, GH2, and GO2. The mixing layers 
among GRP, GH2 and GO2 can be seen from those three contours. The high Oo2 region is the 
flow recirculation region and the flame holder. The pressure contours show a nearly constant 
pressure inside the combustion chamber, which is expected since the jet speeds are subsonic. 
The next three contours show the mass fractions of the three major combustion products; OH, 
H2O and CO2. The high concentration regions represent the flame front. Following is the Cs 
mass fraction contours in which the incipient soot formation starts at the mixing layer. Most of 
the soot in the flame zone and nozzle are burned off by the heterogeneous and homogeneous 
oxidation reactions. Then there is the temperature contours and again the high temperature 
region represents the flame front. The close association of high OH concentration with high 
temperature emphasizes the importance of OH in initiating and sustaining the combustion. 

Figure 4 shows the computed scalar contours for the uni-element LRP/GH2/GO2 thrust 
chamber. The first contours show the LRP jet core and the subsequent trajectories and sizes of the 
droplets stripping from the liquid jet. The rest of the scalar contours resemble those of Fig. 3 
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qualitatively, with some differences caused by the modeling of the liquid jet and droplets. For 
example, the length of GRP jet in Fig. 4 appears to be longer than that of Fig. 3, due to added delay 
of atomization and vaporization processes. The overall temperature of Fig. 4 is also lower than that 
of Fig. 3, caused by the latent heat loss of vaporizing liquid droplets. As such, the computed 
chamber pressure of LRP/GH 2 /GO 2 thruster is slightly lower than that of the GRP/GH 2 /GO 2 
thruster, although the chamber pressure is nearly constant. In addition, the overall computed 
soot concentration appears to be lower in the LRP/GH 2 /GO 2 thruster. 


Table 4 Comparison of chamber pressures. 


Pc, atm 

test data 

34.2 

GRP/ GH 2 /GO 2 model 

35.8 

LRP/ GH 2 /GO 2 model 

33.5 


Table 4 shows the comparison of computed chamber pressures with that of the test. It can be 
seen that the computed chamber pressures from both models compare reasonably well with that 
of the measurement, although it was anticipated that the GRP/GH 2 /GO 2 model predicted 
chamber pressure would be slightly higher. On the other hand, empirical values for such 
parameters as the stripping rate and droplet size distribution were used for the multiphase flow 
model. Since latent heat decreases as chamber pressure increases, it is postulated that states of 
“liquid” and “gas” are indistinguishable under high-pressure rocket chamber environments. This 
hypothesis seems to agree with the observation - a faster rate of gasification is associated with 
elevated chamber pressures - of the combustion time for a RP-1 droplet. As such, the solution 
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from the gaseous phase model is probably as justifiable as that of the multiphase flow model, 
under the rocket operating conditions. 

The discussion herein indicates that the proposed thermophysics model is validated to support 
the design of kerosene-fueled injectors and combustors. It also needs to be pointed that the 
reasonable agreement of the chamber pressures in Table 4 is attributed to the proposed 
thermophysics characterization rather than the choice of single-phase or multiple-phase flow 
options. 

NASA MSFC 40k Kerosene Fueled Thruster 

The proposed kerosene combustion thermophysics model is incorporated with FDNS for 
computing and comparing the predicted radiative heat fluxes with those of a measurement for a 
kerosene fueled thruster plume. The radiometers of the test are located downstream of the 
thruster exit plane and the line-of-sights are perpendicular to the plume centerline. Figure 5 
shows typical computed scalar contours for temperature, CO 2 , H 2 O, and Cs mass fractions. The 
four vertical lines across the plume indicate the line-of-sights of the radiation measurement 
where the first vertical line from the nozzle exit plane is designated as the line-of-sight of 
radiometer Qi, and so on. The line-of-sight of radiometer Q 2 is closest to the computed plume 
Mach disc. The computational model is an axisymmetric formulation and the computational 
domain is five times longer than that is shown to include enough plume for wide-angle radiation 
calculations. Details of injection layout in the combustor such as the multiple fuel and oxidizer 
injection ports, a central fuel jet, and a wall film coolant jet are included in the computation. The 
scalar contours in Fig. 5 show striations of those jets converging near the throat, expanding in 
nozzle and thrusting into atmosphere. The interaction of the exhaust plume with the initially 
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quiescent air creates an eddy-like mixing layer where the remaining reactants afterbums. The 
film coolant essentially forms a protective, sooty layer along the thruster wall and afterbums in 
the outer mixing layer along with other reactants, as shown in the Cs contours. Also shown in 
the Cs contours is a central soot core resulting from the central fuel jet injection. The triple 
shocks (lip shock, Mach disc, and reflect shock) are clearly seen in the temperature contours. 
Overall, the computed thermo-flowfield appears to be reasonable. These scalar contours are 
plotted because CO 2 , H 2 O, and Cs are major radiators of a keosene fueled thmster plume. The 
local concentrations of these radiators and temperature decide the strength of the radiation. 


Table 5. Comparison of computed radiation with that of test measurement. 



0 

4 

180 

4 

5.5 

Q, test 

10-12 

5.4-6.0 

75-95 

>70 

Q, without rxn 

120.7 

2.9 

55.7 

34.0 

Q, irrev. rxn 

8.4 

1.9 

26.5 

25.4 

Q, rev. rxn 

14.1 

5.4 

131.9 

117.3 


Parametric studies are performed to study the effect of homogeneous soot oxidation on the 
radiation signals. Table 5 shows the comparison of computed radiation with that of the test data. 
The hot-firing of a rocket thruster is transient in nature and the radiation signals are taken when 
the chamber pressure reaches an approximate constant. As a result, the plume boundary layer is 
unsteady and the goodness of the radiation measurement often depends on the view angle of the 
radiometer. It can be seen from the view-angle row of Table 5 that only radiometer Q 2 has a 
wide view-angle (180 deg) and is least likely to be influenced by the unsteady motion of the free 
shear layer. On the other hand, the measurements from other three radiometers are of the 
narrow-angle category thereby are more susceptible to the unsteady plume motion. It is also 
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noted from post-test examination that there was some shrinkage of the throat due to soot 
accumulation, resulting in some slight drifting of the chamber pressure. 

From Table 5, the thermophysics model using the reversible homogeneous soot oxidation 
reaction produces the best comparison with the wide-angle data and reasonable comparison with 
those of the narrow-angle measurements (two out of three). The overprediction at Q 3 is 
attributed to the combination of the narrow view-angle and the unsteadiness of the free shear 
layer. On the other hand, when the original irreversible form of the homogeneous soot oxidation 
reaction is used, soot bums off quickly and all the radiations are underpredicted. Without the 
homogeneous soot oxidation, the model predictions disagree with the measurement in all 
respects. The result from Table 5 shows that the proposed thermophysics model provided 
reasonable plume properties. It also emphasized the importance of including the reversible 
homogeneous soot oxidation reaction in the kerosene combustion kinetics model. 

Conclusion 

A thermophysics characterization of kerosene combustion is developed for incorporation with 
advanced computational fluid dynamics methodologies for kerosene-fueled rocket engine design 
calculations. The model consists of a surrogate fuel model representing kerosene/RP- 1 fuels, a 
thermodynamics model anchored with the heat of combustion, and an efficient quasi-global 
kinetics model describing the combustion processes, including soot formation and oxidation 
reactions. The rationale of the modeling approach is given. The benchmark results indicate that 
the proposed thermophysics model can be used to support the kerosene fueled thrust chamber 
design applications and plume definitions. 
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Fig. 3 Computed scalar contours for an unielement GRP/GH 2 /GO 2 injector/combustor: atRp, 0 - 1 ; 
OfcRP. 0-1; 0 Ch2> 0-1; 002 ^ 0 - 1; P, 4.7 - 34.2; ocoh, 0-0.12; 0tH2o, 0-0.87; otco2» 0 - 0.39; Ocs> 0 
-0.17; and T, 251-3528. 



Fig. 4 Computed scalar contours for an unielement LRP/GH 2 /GO 2 injector/combustor: oclrp, 0 - 1 ; 
OcRP, 0 — 1 ; (Xh2» 0 - U cto2> 0 - 1 ; P, 4.7 - 34.2; 0(oh» 0 - 0. 12; (Xh207 0 — 0.87; C(co2> 0 - 0.39; Ocs^ 0 
-0.17; and T, 251 -3528. 


Fig. 5 Computed scalar contours for a kerosene fueled thruster and plume: T, 186 - 6089; Oco 2 . 0 - 
0.5; (Xh 20 ^ 0 - 0-3; and Ocs, 0 - 0.3. 



